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Josephson plasma resonance has been introduced recently as a powerful tool to probe interlayer 
Josephson coupling in different regions of the vortex phase diagram in layered superconductors. In 
the liquid phase, the high temperature expansion with respect to the Josephson coupling connects 
the Josephson plasma frequency with the phase correlation function. This function, in turn, is 
directly related to the pair distribution function of the liquid. We develop a recipe to extract the 
phase and density correlation functions from the dependencies of the plasma resonance frequency 
ojp(B) and the c axis conductivity <r c (B) on the afc-component of the magnetic field at fixed c 
-component. Using Langevin dynamic simulations of two-dimensional vortex arrays we calculate 
density and phase correlation functions at different temperatures. Calculated phase correlations 
describe very well the experimental angular dependence of the plasma resonance field. We also 
demonstrate that in the case of weak damping in the liquid phase, broadening of the JPR line 
is caused mainly by random Josephson coupling arising from the density fluctuations of pancake 
vortices. In this case the JPR line has a universal shape, which is determined only by parameters 
of the superconductors and temperature. 



I. INTRODUCTION 



Recent studies of Josephsen plasma resonance (JPR) absorption in highly anisotropic layered high-temperature and 
organic superconductorstTcj have proved that this tool provides unique information on Josephson coupling of layers 
and on the structure of vortex phases in the presence of an applied magnetic field. JPR in layered superconductors 
was observed by applying a microwave electric field oriented along the c axis which excites charge oscillations between 
layers. In the superconducting state interlayer charge transfer is provided by tunneling of the Cooper pairs. In 
the abseace-of a magnetic field the JPR frequency, oj p , at zero temperature is determined by Josephson interlayer 
couplingErllj: 

^(B = 0)^ 2 = ^Jo = 4^^^ W 
e c $o e c X 2 c ecTKb 

where Jo = c<I>o/(87r 2 sA 2 ) is the Josephson current, e c is the high frequency dielectric constant for electric fields along 
the c axis, X a b and A c are the London penetration lengths for supercurrent in the ab plane and in the c direction 
respectively, 7 = X c /X a b is the anisotropy ratio and s is the interlayer spacing. For the most investigated HTS 
compound B^S^CaC^On (Bi-2212) the JPR frequency v v = lu p /2tt at low temperatures is in the range of several 
hundred gigahertz. At nonzero temperature, in the presence of thermally excited quasiparticles, Eq. p) is valid as 
well, if the JPR frequency is well below the scattering rate of quasiparticles, because under this condition only the 
Cooper pairs contribute to plasma oscillationaljllj. Then, the JPR frequency drops with temperature proportional 
to 1/A^(T). This allows one to observe resonance absopation as a function of temperature at a fixed frequency of 
uniform microwave [electric field oriented along the c axiala. Recently the JPR frequency u>o(T) was also measured by 
sweeping frequencyO and the value wo(0)/27r = 125 GHzj-was obtained for slightly uHderdoped Bi-2212. The same 
JPR frequency was recently obtained for a Bi-2212 whiskertS Transport measurements^ on the mesa fabricated from 
this wisker gave the critical current density Jq(0) — 1200 A/cm 2 . This allows one to extract the dielectric constant 
e c from Eq. (@), e c » 12. 

A DC—magnetic field applied along the c axis, penetrates inside the superconductor in the form of pancake 
vorticeal3~Ej. In an ideal material at low temperatures pancakes would form straight lines and have no influence 
on interlayer coupling. Due to very weak coupling, pinning and thermal-fluctuations easily misalign pancakes in dif- 
ferent layers, which leads to suppression of effective Josephson coupling 2 ^ and, as a consequence, to reduction of the 
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Josephson plasma resonance frequencyciacj. In real materials pancakes in neighboring layers are always misaligned 
due to pinning in the vortex glass phase (below the irreversibility line) and due to thermal disorder in the vortex 
liquid phase. Misaligned pancakes induce a gauge-invariant phase difference <^„, rl+ £(ip)-|between neighboring layers n 
and n + 1 which causes suppression of the JPR frequency according to the relationE3a: 

lj 2 p (B 1 T)=lo 2 (T)C(B,T), C(B,T) = (cos^„,„+i(r)), (2) 

where (. . .) means thermal average and average over disorder caused by pinning centers. 

The JPR frequency is given by the simple expression (^J) and is not sensitive to phase dynamics only if l^(B,T) 
is much higher than the characteristic frequencies u> ps of interlayer phase, slips induced by vortex motioned. This 
frequency can be estimated from the ab component of static resistivity! 2 ^, uj ps ~ (2 -j- 4) • 10 8 [l/s]-T[K]-p a h. We 
will show that, more directly, uj ps can be extracted from both the c axis conductivity and plasma frequency. These 
estimates show that typically oj ps in the liquid region varies from 10~ a uj p at T w 40 K to 10~ 2 cl> p at T k, 70 K. Thus, 
the plasma frequency is determined by a snapshot of vortex positions. For this reason JPR study provides information 
on the equilibrium properties of the vortex lattice though it is dynamic in nature. 

The plasma frequency is mainly sensitive to the average value of cos (r) even if its local variations are very 

large, provided the random part of cos</? n! „+i(r) changes rapidly in space. In this case these rapid variations are 
averaged out by the smooth external oscillating electric field. Space variations of cos</3 nin+ i(r) result in mixing 
of phase collective modes with different momenta and frequency. Such a mixing is responsible for resonance line 
broadeningc3'Ej. Mixing is determined by Josephson coupling and is proportional to / 2 /A 2 , where Aj = 7s is the 
Josephson length, and l v is the characteristic length of the phase correlation function 

S(r, T, B) = (cos[ ¥ >„, n+1 (r) - <^ n+1 (0)]) . (3) 

Averaging in space is effective for l v <C Aj. This condition is fulfilled for fields B 3> Bj = $o/A 2 7 in a vortex phase 
with strong disorder along the c axis, when l v w a. Such a situation is realized in the vortex liquid phase where 
interlayer correlations of vortex positions are very weak and also in the Bean critical state obtained by sweeping 
magnetic field after zero field cooling. l v is of the order of a also in the vortex glass phase obtained by cooling in 
high magnetic fields above the "second, peak" field (ss 300 — fiOfl G), where the three-dimensional vortex lattice is 
strongly disordered according to /i + SRE3 and neutron scatteringclll measurements. In this decoupled phase, interaction 
of pancake vortices in different layers is very weak, see discussion in Ref. |3l| . 

Therefore, the key parameter which determines field and temperature dependence of the JPR frequency is C(B, T). 
This allows the use of JPR as a sensitive probe of the vortex state. The field dependence of uj p allows one to observe 
JPR absorption as a function of magnetic field by sweeping the magnetic field B at fixed frequency u> of applied 
microwave electric field. In this way the resonance magnetic field B r was obtained as a function of temperature and 
frequency lu, ui p (B r , T) — u>, see, e. g., Refs. Tsui et ala found the field and temperature dependence of uj p in 

the vortex state of Bi-2212 by measuring B r (ui, T) at several frequencies and using field sweeping (FS) after zero field 
cooling. They found that the resonance frequency depends nonmonotonically on temperature at fixed field and has 
a maximum at the irreversibility line. Measuring B r (uj,T) in the field cooling (FC) and FS modes, Matsuda et al.a 
demonstrated that the plasma frequency is history dependent below the irreversibility line. Recenly the Josephson 
plasma resonance has also been measured in the vortex crystal state at small fields£3Cj. Field dependence of the JPR 
frequency in this phase is well described by the model based on the suppression of the Josephson coupling due to 
elastic fluctuations of pancake vorticescil. 

In the vortex liquid phase, for the magnetic ficlji^pBp oriented along the c axis, the field and temperature dependence 
of JPR frequency is approximately described ast 



(B Z ,T)^^(T)B- V T~ 1 (4) 



with v being slightly smaller than unity (« 0.8). 

The influence pf vortices on interlayer coupling in the vortex liquid was also probed by direct c-axis transport 
measurements£3~E3. In particular, pancake misalignment suppresses the Josephson current J C (T,B Z ) as compared to 
its zero field value J (T), as J C (T,B Z ) « J (T)C(T, B z ), so that J C (T,B Z ) is proportional to c^(B 2 ,r). The field 
dependence J C (T, B z ) cx Bj 1 , similar to Eq. ( ||) was indeed found for the c-axis critical current) 34 ' 



Behavior of JPR in the liquid phase described by Eq. (|4|) was explained in Ref. 25 in the framework of the high 
temperature expansion with respect to Josephson coupling. This approach provides a general relation, connecting the 
average cosine factor C(B,T) with the phase correlation function S(r, T, B) defined by Eq. (|3|), 



*JWT,B.)-,™. (5) 
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where thermal average (. . .) in Eq. (||) is now done without the Josephson coupling, Ej — Eq/X 2 is the energy of 
Josephson interlayer coupling per unit area, Eq = ^>Qs/16n 3 X 2 b , and 

f=^-JdvS{v,T,B z ) (6) 

is a dimensionless universal function of order unity. The relation (|J) was derived assuming that: a) (cos[i/? njTl +i(r)]) = 
if the Josephson coupling is absent, and b) S(r,T,B z ) is determined by vortices only and thus its correlation length 
l v ~ a. Such a perturbative approach is valid if the Josephson energy in the area I 2 , is small compared with T, i.e. 
Eolu>/Xj ~ E a 2 /Aj <C T. This condition is fulfilled in magnetic fields B z » Bj(E /T) and, under such condition, 
one gets uj p (B,T) <C luq(T). For Bi-2212 crystals with 7 of the order of several hundred the high temperature 
expansion approach is valid practically at all fields above w 100 G. 

The same physical picture provides prediction for dependence of the JPR frequency on the in-plane field B x . In the 
decoupled phase at B z 3> Bj the in-plane field produces an additional linearly growing term in the phase difference, 
2irB x ys/$o, which accounts for the phase difference due to the vector potential corresponding to B. x . Therefore, the 
dependence of w 2 on B x is simply determined by the Fourier transform of the phase correlation function S(r,T, B z ). 
It was also assumed in Ref. |2j] that S(r,T,B z ) can be approximated by a Gaussian, S(r) = exp(— nB z r 2 In 
this case the dependence of JPR on the in-plane field can be presented in the analytical form: 

u 2 p {B x ,B z ,T) « u$(B„T)exp(-fin?Bl/$ B,). (7) 

The exponential factor is the Fourier component of ^(riri Such dependence was found to describe fairly well the 
experimental angular dependencies of the resonance fieldcm. 

In this paper we develop a detailed theory of the Josephson plasma resonance in the liquid phase based on perturba- 
tion theory with respect to the Josephson coupling. This approach is justified because ui 2 /uj 2 = C(B, T) is well below 
unity in this state. We will extend the analysis of Ref. ^5] by taking into account spin- wave type phase fluctuations in 
addition to the effect of vortices on JPR. In this approach S(r,T,B) — S r (r,T)S v (r,T, B), where S r (r,T) is deter- 
mined by spin- wave type phase fluctuations, while S v (t,T, B) is the phase correlation function induced by vortices. 
We will calculate S r (r, T) and express S v (r, T, B) via the correlation function of vortex density in the system without 
Josephson coupling. We show that, at temperatures above the melting point, the coupling between two dimensional 
pancake liquids in different layers is very weak and can be neglected. If, in addition, pinning inside the layers is 
neglected, then the density correlation function coincides with the density correlation function of the two-dimensional 
one-component Coulomb plasmaEJ. Accurate analysis shows that the Gaussian dependence of S v (r) used in Refs. |25| 
and |] actually exists only at r < a; at large distances S v (r) decays exponentially, S v (r) oc exp(— r/l). This means that 
the dependence of JPR frequency on the in-plane field B x is also not exactly Gaussian, and is given by the Fourier 
transform of the phase correlation function, f s (b) ~ (1/a 2 ) / dr exp(ibx / a) S (r) , 

lo 2 (B x ,B z ) = c^(Q, B z )f 3 (B x /y/B z H Q ), (8) 

with Hq = $o/(27rs) 2 . As the function f s (b) is directly connected to the phase and density correlation functions, 
this equation allows us to obtain information on the structure of the vortex state from experimental study of JPR. 
We calculated numerically the two-dimensional density correlation functions using Langevin dynamics simulations 
and used these functions to calculate the phase correlation functions S v (r,T,B) at different temperatures. Influence 
of weak interlayer coupling and pinning is investigated perturbatively. We found that the function f s (b) in Eq. (@) 
calculated from numerical data gives a better description of the angular dependence-pf the resonance field than the 
naive Gaussian ansatz (Q). Part of these results has been published in a short paperca. 

We will investigate also the shape of the resonance line. We demonstrate that in the case of weak damping in 
the liquid phase the plasma resonance line acquires a universal shape due to fluctuation-induced inhomogeneities of 
the interlayer phase difference. This mechanism gives a natural explanation for experimentally observed asymmetric 
broadening of the line. The tail at high frequencies arises due to mixing of propagating plasma modes by random 
Josephson coupling. The low frequency tail is caused by plasma modes localized in the rare large areas where the 
Josephson coupling is suppressed due to fluctuations, similar to the well known Lifshitz tail in the electron density 
of states in disordered semiconductors. These results were published in the short communication |28|. In this paper 
we present their detailed derivations. At high frequences (fields) we calculate the shape of the line using perturbative 
expansion with respect to inhomogeneous Josephson coupling. At frequencies below the resonance frequency we use 
the method of optimal fluctuation to calculate the exponential tail of the JPR line. 

The paper is organized as follows. In Section II we present general thermodynamic relations for a layered super- 
conductor in the mixed state and discuss phase dynamics. Section III is devoted to Josephson coupling in the liquid 
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state. We obtain the high temperature expansion for the parameter C(B,T), derive general relations connecting the 
JPR frequency with the pair distribution functions, and calculate the pair distribution functions numerically. We also 
establish a relation between the plasma frequency and the static c-axis conductivity. In Section IV we calculate the 
shape of the JPR line in the liquid state due to the inhomogeneous broadening mechanism. 



II. GENERAL FORMALISM 



A. Thermodynamics 



We consider a layered superconductor in the vortex state. In the most part of the field-temperature phase diagram 
below T c the relevant degrees of freedom are the vortex coordinates r nv and phase spin- wave type degrees of freedom 
Pnli+iW- The energy functional in terms of these degrees of freedom can be written asEZI: 



(9) 



Here T em {v nv ) is the functional jsdiich accounts for the two-dimensional energy of pancakes and their electromagnetic 
interaction in different layerslI3~EII: 



3~ sm^ni/) 



Eo 

47T 



dkdq 



exp^k ■ (r ni/ - r mv >) + iq(n - m)] 



(10) 



77, 777, 1/, T'' 

where Q 2 = 2(1 — cosq)/s 2 . The contribution T V i n (? n v) accounts for pinning: 



T V mk?nv) = ^2 J drV pin (r)5{r - r nv ), 



(11) 



where V^„(r) is the pinning potential for vortices, see, e.g., Ref. |4^. The last term is the energy of intralayer currents 
associated with spin- wave type of excitations and the Josephson energy. It depends on the regular part of the phase 
difference and vortex coordinates: 



(12) 



I m.n. J 7). 



Here V = (d/dx, d/dy). The mutual inductance of layers, L nm , is: 

\n — m\s 



The phase i/?^ +1 (r) is the singular part of the phase difference induced by vortices at positions r„„ and r n +i v when 

(14) 



— L{q) cos(n - m)q « — exp 
^ 2tt Zs 



A a 6 



L(q) = 



1 



2(l-co S(Z ) + s 2 /A^ 



(13) 



Josephson coupling is absent (1/Aj = 0): 

-i( r ) = X^( r ~ tw) -<j> v (r- r n+ i jV )], 



(v) 

^77,77+1 



where 4> v {r) = arctan(y/x) is the polar angle of the point r. 

The average cosine of the total phase difference is given by the expression: 

T d 

C = l+ — (lnZ) dis , 

NSdEj K 1 ' 

Z = [ D^l +1 (v)Dv nv exp\-^{v nu ,^ n+1 (v)}/T 



(15) 
(16) 



Here 



is an average over disorder (pinning), N is the total number of layers, and S is the total area of the layer. 
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B. Phase dynamics 



Time variations of the phase difference leading to excitations of quasiparticles and the deviations of the quasiparticlc 
distribution function from equilibrium should be accounted for in the time-dependent equation for the phase difference. 
General equations describing phase dynamics are still under discussion and they depend on the mechanism of pairing, 
see Refs. |f3|-|45|. At the phenomenological level, the dynamic equations for the phase difference </5 n ,n+i (r, t) can 
be derived as follows. Introducing the in-plane superfluid momentum p„ = V^> n — (27r/$o)A„, we can express the 
in-plane magnetic field as 

$o 

H_l = - — n z x (X7tp n +1 - sV z p„) 

Z7TS 

where n z is the unit vector in z-direction, V 2 p„ = (p n +i — Pn)/ S - Using this representation the z-component of the 
Maxwell equation cur IB = (47r/c)j + (l/cjffD/dt can be written as 

47T . s c dE z $ /-,9 —„ s f*rt 

c c ot 2ns v ' 

where the current j z consists of the Josephson and quasiparticle contributions, j z — Jo sin tp n ,n+i + jqz- From the 
in-plane component of the Maxwell equation we obtain the independent equation for Vp„ 

(-V2 + a- 6 2 ) v P „ = -gv (j„ + g^) - iv, W-x,n, (18) 

Here e a b is the high frequency dielectric constant for the in-plane electric field. We split the in-plane current into the 
quasiparticle and superconducting components, j = j s + j q with j s = (c < i>o/87r 2 A^ b )p n . In the following we consider 
the case when the in-plane equilibration occurs at time scales much shorter than the inverse plasma frequency. This 
allows us to neglect the in-plane quasiparticle current and electric field in Eq. (|l^). Solving formally Eq. dig) with 
respect to Vp„ and substituting it into Eq. (fl7|) we obtain 

Sin - Aj y^£„ m VVm,m+l + ^ + T-V-Ir = ( 19 ) 

^ Jo 47rJ ot 

To obtain a closed set of equations for (p n . n +i, we have to connect E z with dip n>n +i/dt. In a common Josephson 
junction theae-quantities are connected by the famous Josephson relation d<j) n ,n+i/dt = (2ncs/(&o)E z . However, it 
was shownEra that in atomically layered superconductors this simple relation is violated due to the charging of 
layers. The charging effect leads to the variations of the chemical potential \i n in the layer n, and, in general, the 
time derivative of f n ,n+i can be written as 

—d^-^\ Ez - Vz -)- (20) 

where fi n is related with the variation of the charge p n as 

47reA 2 n , . 

M„ = -Pn (21) 

e c s 

where A^ is the Debye screening length. The charge density is related to electric field by the Poisson equation 

An 

V 2 £ 2 = —p n (22) 

Sc 

Combining Eqs.([20|), (pl|), and ( p2] ) we obtain relation between E z and dip n ,n+i/dt 

dLpn > n+1 - 2nCS (l-\lVl)E z . (23) 



dt $ 

Finally, the quasiparticle current is determined by the quasiparticle conductivity er q2 



. _ $QCT qz dip 7hn+ i 

JqZ ~ 2nCS Ot 1 > 
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The interlayer quasiparticle conductivity in superconductors with d-wave pairing remains nonzero even at T 
and it was found to be w 1 — 2 (kohm cm) -1 in Bi-2212 crystals at 4 K and B = OcJ depending weakly on Be 



Eqs. (|19|), (p3[), and ( 24J ) form a closed set of equations for (p n , n +i- They describe phase dynamics when deviations 
of the quasiparticle distribution function from that at equilibrium can be neglected. Deviations from equilibrium 
and charging effects may lead to the broadening of JPR resonance in addition to the effect of dissipation caused by 
quasiparticles. It is unclear at that moment how important-are noncquilibrium effects because until now they were 
studied only in the framework of simple microscopic modelsc3, which may be unadequite for HTS. 

As we are mainly interested in the effect of pancake vortices on JPR in the follco&apg, we will use simplified equations 
in which we do not account for quasiparticle conductivity^ and charging effectsLL3c3. This simplest form of equation 
for tfn n+iir, t) which describes phase variations at fixed vortex positions is: 



9. 



S+i-^E^Cm+i + sm <;+i + <l+i(r) =^j-Qi V S), (25) 



1 d 



where T) z (t) is the homogeneous in space alternating external electric field, and we assume that tp^\ +1 (r) is time 
independent, because characteristic frequencies of vortex motion are well below the plasma frequency as discussed 
above. For small amplitude variations we expand the phase difference as: 

<^„, n+ l (r,t) = ¥>i°,n+i( r ) +^,„+i(r,t). (26) 

Here <fn\i+i( r ) i s ^ ne solution of the stationary equation. We obtain for the Fourier component of the phase difference 
induced by the external electric field, ip' n n+1 (r,tj), the equation 

2 

^2<n+l + A} ^C+l " COf #i 0) «+l( r )] = ~T^T V ^ t )^ ( 2? ) 

UJ a 47TJo 

u m 

The homogeneous part of the averaged electric field between layers n and n + 1 is 



(£*M) = |^« n+1 ) (28) 



and JPR absorption is determined by the imaginary part of the inverse dielectric function, l/e(uj) = (E z (uj)) /T> z (ui). 
Expanding ip' n n+1 (r) in terms of the eigenfunctions ^ an {v) of the Schrodinger-like equation 

2 

- V L nm V 2 V am + cos[^ 0) n+1 (r)]^ Qn = ^ an , (29) 

LOr. 



we obtain a useful expression for Im(l/e(o;) 

1 7TUJ 

Im 



e{uo) 2e c 



j dr (J2 ^ao(0)Kn(r)S(uJ -u a )Y (30) 



where (...) is the average over configurations of cosy>„°„ +1 (r). This expression will be used in Sec. V to calculate the 
shape of the JPR line. 

III. JOSEPHSON COUPLING IN THE LIQUID STATE 

A. High temperature expansion 

In the vortex liquid phase we can use perturbation theory with respect to Josephson coupling, which relates the 
cosine of the total phase difference with the phase correlation function S(r) as in Eq. (||). Splitting the total phase 
difference <p n! „+i(r) into vortex, fn n +i( r )i and regular Lp^\ l+1 (r) contributions, we obtain: 

( C0S [<Pn,l+i + fnjn+il) = Tpf / dr{cos[tp ntn+1 (r) - <p n ,n+i(0)]) = drS v {r)S r (r), (31) 



G 



where the correlation function of phases S v (r) = (cos[<^^ +1 (r) — fnii+ii^)]) ^ s determined by vortex positions, 
S v (r) = (Z- 1 J Dv nv cos[^ +1 (r) - <^J +1 (0)] exp{-(l/T)[^ em + F pin ]}) dis , 
Z v = j Dr nv exp{-(l/T)[F em + F pin }}, 

and Sy(r) = (cos[<£^ +1 (r) — Pn n+ii®)}) describes the effect of spin- wave type phase fluctuations (phase difference 
fluctuations in the absence of vortices), 

Sr(r) = Z- 1 I D^ n+1 cos[^ +1 (r) - <p% +1 (0)] exp[-(^/2T) $>„ m V<^ +1 • V^ ro+1 ], 

nm 

Zs= I C^ +1 exp[-(So/2r)^L„ m V^ +1 -V^ m+1 ]. 

nm 

In the lowest order in Ej, spin- wave type of fluctuations and vortex coordinates degrees of freedom do not interact. 
In the Gaussian approximation for spin-wave type of phase fluctuations, we obtain: 

T /■ I — cos(kr) \ (U\ 2a T 



SAV) = ° XP {- W J dMq im^ ) = {r) - ° =2^- (32) 

We cut off the logarithmic divergence in the integration over k is at k w l/£ a &, where is the superconducting 
correlation length. As we show below, a particular value of this cutoff has no influence on observable quantities 

B. Relation between the phase correlation function S v (r) and the density correlation function of the liquid 

Next we express the phase correlation function S v (r,T,B) via the density correlation function of pancake vortices 
induced by the B z component of the magnetic field. Let us consider the phase difference between points r and in 
the layer n induced by pancakes with density p„(R) = Yl u ^(R — r n»)- Denoting the phase in the layer n at the point 
r by (j> n {v) we obtain 

0„(r)-0„(O) = yrfRp„(R)/3(r,R), (33) 

where 

/?(r,R) = ^(r/2-R)-^(-r/2-R) (34) 

is the angle at which the segment connecting points — r/2 and r/2 is seen from the point R = (X, Y). (3(r, R) has a 
jump of 2tt when point R intersects segment [—r/2, r/2]. Choosing the x axis along the segment [—r/2, r/2], we can 
represent /3(r, R) as 

Yr 

/3(r ' R) = arCSin [(^ + rV4)^-(^]V 2 ' (35) 

For the phase difference <pi"i+i( r ) — Vnn+iW between layers n and n + 1 in the presence of both components B x 
and B z we obtain 

r{n+l)s 

^ + i(r)-^ + i(0) = ^(r) + (27r/$ ) / dzL4 z (r,z)-^(0,z)], (36) 

*„(r) = J rfR[p„(R)-p„+i(R)]/J(r,R). (37) 

Here ^4 z (y) = —B x y is the vector potential for the -Bj, component of magnetic field. We assume here that B x is large 
enough so that it penetrates almost freely into the sample and therefore has no influence on the phase fluctuations. 

In the following we will use the Gaussian approximation for the random function $„(r), i.e. we neglect irreducible 
higher-order correlation functions. Then 
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(cos[*„(r)]) « exp[-F„(r)], F v (r) = ([<P v (v)}) 2 /2. (38) 

Therefore the vortex-induced phase fluctuations are characterized by the the correlation function 

S v (r) = exp[-F„(r)] co$(2TrsB x y/$ ). (39) 

Using Eqs. ( |37|) and (|3|) we connect F v (r) with the following density correlation function 

K(v) = (l/2)([p„(r) - Pn+1 (r)][p n (0) - p n+1 (0)}) 

= ^o(r)-^i(r) (40) 

by the relation 

F v (r) = J dRdRi/v (R - Ri)/3(r, R),3(r, Ri), (41) 

where K n (r) = (6po(v)Sp n (0)} and Sp n (r) = p n (r) — n v is the deviation of vortex density from its average value 
n v = B z /$ in the layer n at point r. We will use the representation 

JTW-mW+nWr), -./«*■> ~1. (42) 

In the case of uncorrelated pancake arrangements in neighboring layers h(r) is the pair distribution function of the 
2D liquid. Using identity J drK(r) =0we rewrite Eq. (41) as 

F v (r) = -(n 2 v /2) J dndR/i(n)|/J(r, R - n/2) - /?(r, R + rx/2)] 2 (43) 

Performing the integration over R and averaging over orientation of i"i we finally obtain 

F v {r) = -irn 2 v J dnrxhin) J(r, n), (44) 

where the universal function J(r, ri) is given by 

J(r, n) = 2?rr 2 [ln(r x /r) + 2] , n > r, (45) 
J(r, ri) = 47r [2rri - r\ - (r 2 /2) ln(r/ri)] , t\<t. 

The function J(r, ri) and its first three derivatives are continuous functions of r and r\ but the fourth derivative has 
jump 87r/r 2 at r = r%. Eqs. ( ^i[ ) and (45) represent a very important general relation connecting phase fluctuations 
with the vortex density fluctuations in the Gaussian approximation for the pancake liquid. The function F v (r) 
increases linearly at large r, 

poo 

F v (r) mr/lo, a/l a = -87r 2 / dxx 2 h(x), r > a, (46) 

Jo 

indicating exponential decay of phase fluctuations at large distances. Here we introduced the dimensionless distance 
x = r/a with a — n v . For small r the function F v (r) increases approximately quadratically with r: 

F v (r) w 7T (r/a) 2 (in- + 2 + C 2 J , C 2 = -2ir J dxxh(x)lnx, r < a. (47) 

Note that depends on the temperature and pinning. The temperature dependence is determined by the dimen- 
sionless parameter T/Eq if the effect of pinning is neglected. Respectively, the function F v (r/a) is also temperature 
and pinning dependent. 
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C. Effective Josephson coupling and plasma frequency 



The problem of calculating of C(B,T), which determines the Josephson plasma frequency is now reduced to cal- 
culation of the integral ( f4~4| ) with pair distribution functions of the, liquid, h(r). These functions are not available 
analytically. They can be calculated from Monte Carlo simulationstiH or (approximately) using the density functional 
theorycj. As the first approximation we consider uncorrelated two dimensional liquids in layers. For weakly coupled 
superconductors the density correlations between the layers and weak pinning can be taken into account perturba- 
tively. In the 2D approximation the vortex system-is equivalent to the one-component two-dimensional Coulomb 
plasma, which was studied extensively in the pastfiH. The pair distribution function /i2u(r, B, T, Eq) within this 
approximation has an exact scaling property: 

h 2D (r,B,T,E ) = h 2D {r/a,T/E ) (48) 

i.e., it depends on magnetic field only through the spatial scale. As follows from Eq. ([44| ) this scaling property is 
also transferred to the function F v (r), F v (r) — F v ( r/a, T /Eq). If, in addition, fluctuations of the regular phase 
are neglected then at B x = we obtain from Eqs. ( pl| , p9| ) the expression (|^) for C. Correlation properties of the 
liquid determine the field independent universal function J{T/Eo) defined by Eq. (^). Therefore, the scaling property 
C oc 1 / B z would be exact if we neglect (i) coupling between pancake liquids in different layers (ii)pinning inside the 
layers, and (iii) regular phase fluctuations. All these effects are typically small in the pancake liquid phase. 

Now, we further develop this approach and include fluctuations of the regular phase, magnetic interlayer coupling 
and pinning. First, we include the effect of regular phase fluctuations. As the temperature approaches T c the role of 
these fluctuations of regular phase progressively increases. The correlation function of regular phase S r (r) is given 
by Eq. (|||). We substitute ( p^ ) into Eq. (^) and take into account that the Josephson coupling, Ej oc A^ 2 , and the 
plasma frequency at zero field are renormalized by the phase fluctuations. Their suppression is determined by the 
cosine factor at B = 0, C(T), which was estimated in Ref. ^ as 

C(T) = {UI\jY- (49) 
Using the scaling property of S v (r) we obtain for magnetic field along the c axis: 



where 



l-a 



cos[<4 +1 + = /.(T)-^ ) . (50) 



f s (T) = 2ir dxx 1 - 2 " 
Jo 



S v (x,T) 



again depends only on reduced temperature T/Eq within the 2D approximation. Therefore, regular phase fluctuations 
reduce the power index in the fi eld dependence of C and make it temperature dependent. Numerical calculations, 



described below in Section |III Ej , give f s (T) « 1.0 - X.75T/(2irE ) at a = T/(2ttE ) < 0.1. Using Eqs. fl3l|), (|32|) and 
( [i"l| ) we obtain the JPR frequency for arbitrarily oriented magnetic field: 



^( B „ B „r, = „ S( r)§(|) , -" / .(^,r), m 

poo 

f s {b,T) = 2n dxx^ 2a S v (x,T)J {bx), (52) 



where Jq(x) is the Bessel function. 



D. Influence of interlayer coupling and pinning. Perturbative estimate. 



We consider now the influence of weak interlayer magnetic coupling and weak pinning. In both cases it is necessary 
to calculate a correction to the density correlation function and substitute it into Eq. (44). The starting point of the 
calculation is the thermodynamic expression for the density correlation function 



K n -m.{r - r') ee (5p n (r)Sp m (r')) 
1 



Dr n „6p n (r)Sp m (r') exp 



J~2D + ^roc 



(53) 
(54) 
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Here we split the electromagnetic part of the energy (|10|) into the intralayer two-dimensional part (T 2 d) and the 
interlayer magnetic coupling (F mc ) part. It is convenient to represent all energy contributions in continuous form: 

T 2D = nE o y f dvdv 1 In —2—5p n {r)8 Pn {r') (55) 
„ J |r-r'| 

fmc = 2 X] / drdr ' v ™d r - r', n - m)5p n (r)Sp m (r') (56) 

^ / drF p in(r)5p„(r) (57) 



pin 

with 

(Vpin(r)V p in(0)) = W(r). (58) 

Magnetic interaction between pancake vortices, 14i C (r, ra), has been considered in Rcfs. |l9| ^l[ At small distances 
r, ns <C X a b a simple analytical expression can be obtained for n > 



U mc (r,n) = 



\J r 2 + (ns) 2 — ns — ns In 



•y/r 2 + (ns) 2 + ris 



2ns 



(59) 



If both f mc and J-pi n are neglected then K n - m (r) — K 2 D(v)8 nm , with l<r 2 £>(r) being the density correlation function 
of a 2D Coulomb plasma, K 2 r>{v) = n v S(r) + n 2 /i2zj(r). 

First, we obtain the correction to the density correlation function induced by magnetic coupling, A mc _fir„_ m (r — r'). 
Expanding ( 53) with respect to T mc we obtain: 

A mc lf n _ m (r) = -~ f dridr 2 V mc (r + r 1 -r 2 ,n-m)K 2D (r 1 )K 2D (r 2 ). (60) 

Interplanc phase correlations are determined by the density correlation function in neighboring layers, \m — n\ = 1 
(see Eq. (flo| ) ). The integral ( |60| ) decreases rapidly at r ^S> a, because the main contribution comes from the distances 
ri,r 2 of the order a and J dvK 2 u(Y) = 0. For s <C r < A a &, the interaction potential for pancakes in neighboring 
layers V mc (r, 1) increases linearly with r, U mc (r, 1) w (7rsi?o/2A 2 b )r and the correction to Kx(r) is: 

A mc ^!(r) = -£^= / drxrfr 2 (|r + r x - r 2 | - r) ^(ri^fa). (61) 

^ab 1 J 

Using the scaling property of K 2 d ( r ) i we get a simple dimensional estimate of this integral 

A^K^/nl = {saE Q /\l b T)F mc {r/a), (62) 

where F mc (x) is a dimensionless function of order unity. An exact shape of this function can be obtained by numerical 
evaluation of the integral in Eq. ( |6l| ) . Using the scaling estimate ([32|) we obtain the correction to C given by Eq. ( p30| ) 
due to magnetic coupling, A mc C: 

"roc^ - Jmc ^ y2£3/2 

with / mc ~ 1. This correction is small because of the small parameter s/X a b ~ 0.01 in atomically layered supercon- 
ductors. 

Let us turn now to the pinning-induced correction to the density correlation function, A p i n i^„_ m (r — r') = 
Api n -ftTo( r — r')5 nm . Expanding ( |5^ ) with respect to .F p i n we obtain in the second order with respect to pinning 
potential: 

A pin X (r-r') = (63) 
^2 / d ri dr 2 W (n - r 2 ) [(<5p(r)<yp(r , )«p(ri)*p(r 2 )) - K 2D (r - r')X 2D (r a - r 2 )] . 

This correction is determined by the quartic correlation function of density. Such a perturbation approach is valid for 
weak pinning when the pinning energy of a single vortex, E p s» y/W(0), is small in comparison with the temperature. 
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An important observation following from Eq. (|6^) is that, in this temperature interval, the first pinning correction 
exactly coincides with the correction induced by the extra interaction potential W (ri — r^) fT. The range of this 
potential coincides with the correlation length of pinning, r p , that is much smaller than the intervortex spacing. 
It means that this potential has almost no influence on the density correlation function at E p -C T, because the 
contribution to the correction comes only from the rare configurations when two particles approach at distance of 
the order of r p . We can conclude, therefore, that influence of pinning on the pair distribution function of the liquid 
is negligible at T > E p . This conclusion is in agreement with the more detailed study of influence of point pinning 
on correlation properties of the pancake liquid using the replica approachEil. Note, that strong pinning sites with 
E p « Eq , such as columnar defects produced by the heavy ion irradiation, can not be treated perturbatively and may 
strongly influence density correlations of the pancake liquid. 



E. Phase correlations in pancake liquid regime. Numerical calculations 

In the decoupled liquid phase correlations can be calculated using the pair distribution function h(r) of the two- 
dimensional vortex system. The phase correlation function F v (r) is connected with h(r) by relation (|44|). To study 
the behavior of phase correlations in the liquid state we performed numerical simulations of the motion of N v = 900 
point vortices by numerical solution of the set of Langevin equations for the vortex coordinates Ri: 

^ = f v (B< -Rj) + hi (t). (64) 

Here £ u (r) is the intervortex interaction force, fuif) is the Langevin force, (fLia{^)fLjp{t)) = 2T6 a p5ij5(t). The 
quantity 2ttEq — 2s<f>g/(47rA a h) 2 , and the lattice constant ao = y2/ \/3n v of the ideal lattice are taken as the units of 
energy and length. To facilitate numerical calculations we model the real infinite-range 1/r interaction force by the 
finite-range interaction force f v (r) = (1/r) (l — r 2 /r 2 ut ) , cut off at the finite length r cut . To minimize the effect of 
cutting we have chosen the cutting length r cut — 4.33 to reproduce the shear modulus 6*66 = B$o/ (8ir\ a i,) 2 for the 
uncut interactions (i.e., for r cut — oo). We solve the Langevin equations (]64|) at different values of T ranging from 
0.006 to 0.1. Simulations give a melting transition at T m w 0.007, in agreement with the value of the melting point for 
the two-dimensional Coulomb plasmaE3. A set of pair distribution functions, h(r), has been generated using averaging 
over a large number ( w 2- 10 4 ) of vortex configurations (see Fig. Q). For each h(r) we calculated the phase correlation 
function, S(r), integrating numerically Eq. (f44|). Several phase correlation functions are presented in Fig. |^. Note 
that the phase correlation function has surprisingly weak temperature dependence. To facilitate comparison with 
experiment, we also calculated the function f s (b,T) that determines the angular dependence of the plasma frequency 
via Eqs. (|5l|) and ([52|). This function is plotted in Fig. [|. We also calculated typical lengths, characterizing decay of 
S(r): lo defined by Eq. (Q), R, and R%, where the last two lengths are defined by relations 

/•oo />oo 

R 2 = rdrS{r), R\ = / r 3 drS(r)/R 2 . (65) 
Jo Jo 

Temperature dependence of the phase correlation length Iq is shown in Fig. |^. For comparison we also show two 
points extracted directly from the phase correlation functions, which were generated using Monte Carlo simulations of 
the uniformly frustrated XY modeE3. We see that the results of the two models are in perfect agreement. We found 
that, in spite of considerable change in shape of h(r) with increasing temperature, the phase correlation function does 
not change much. In particular, we find that Iq drops from 0.3 ao to 0.25 ao (see Fig. |^), R 2 w (0.12 — 0.13) a 2 , and 
R\ w (0.49 — 0.55) Oq, (Ri ~ 2R ) in the investigated temperature range. 



F. Extracting the phase correlation and pair distribution functions from angular dependence of the plasma 

frequency 

Measurements of dependence of the resonance frequency uj p on B x at fixed B z and T allow one, in principle, 
to obtain the phase correlation function S(r) and the pair distribution function h(x), i.e., to extract quantitative 
information about correlations in the vortex liquid. Using Eq. ( |5l"| ) the function f(b) — f s (b)/f s (0) can be found from 
the dependence lo p (B x ) as: 

u 2 p (B x = bVW;/2n S ) 

m = o) • (66) 
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Usually the angular dependence of the resonance field B r (9) is measured at fixed microwave frequenc; 
u p (B x ). This dependence also can be used to extract f(b): 



f 



b = 



27TsB rx (6) ' 




B rz (9) 


l-a 


^$ Q B rz {6)_ 




B rz (0)_ 





rather than 



(67) 



The function F s (x) may be found by the reverse Fourier transform from Eq. ( |52| ) , 

-f s (0,T)x 2a r°° 



F s (x,T) 



In 



dbf{b)bJ (bx) 



2vr Jo 

In the final stage the density correlation function h{x) according to Eq. ( E3) is given by the relation: 



h(x) 



1 



87TX 3 dx 



dx 



d 3 F s (a 
dx 3 



(68) 



(69) 



Unfortunately it includes high order derivatives which makes the procedure very difficult for practical realization. 

It is more practical to compare the Fourier transform of the phase correlation function fib), which can be directly 
extracted from the angular dependence of resonance field via Eq. (|67|), with the theoretically predicted function 
obtained by numerical simulations (Fig. |^). To make this comparison we used angular dependence of the resonance 
field for an optimally doped Bi-2212 sample at 40 K from Fig. 3 of Ref. [§ We plot in Fig. | the function B rz (9) /B rz (0) 
(a at this temperature can be neglected) vs 2irsB rx {9)/ ^<& G B rz (9) and the theoretical f(b) at T/2nE Q = 0.032, 
roughly corresponding to 40 K. One can see that the experimental f(b) is very close to the theoretical one. This 
indicates that phase fluctuations in real Bi-2212 crystals are well described by the pancake liquid model. 



G. Relation between the plasma frequency and the c-axis conductivity in vortex liquid 

In this part we will establish a relation between uj p and the c -axis conductivity a c in the vortex liquid which allows 
us to estimate the characteristic frequency of phase slips, lu ps , from experimental data for uj p and er c . For the plasma 
frequency we obtain: 

w 2 p {B) = (4^ J$/e c T) J dvS{v, B, t = 0), (70) 

S(r,B,t) = ( (r,*)-v> n> n + i(0,0)]) ) (71) 

The c-axis conductivity below T c consists of the quasiparticle and Joscphson contributions. The conductivity due to 
the interJayer Josephson currents is finite in the liquid phase because of the phase slips induced by pancake thermal 
motioncj. When the temperature is decreased at fixed field, the system smoothly crosses over from the region of 
dominating quasiparticle transport to the region of dominating Josephson transport. We consider the latter region, 
where the conductivity is related to the correlation function of the Josephson currents via the Kubo formula 



cj c = (sJ 7T) J dtj dr(sin(p n ,„ + i(0,0)sin^ n ,„ + i(r,i)) = {sJ£/2T) J drdtS(r,t). (72) 
Now the characteristic frequency uj ps involved in the c axis conductivity may be estimated using the relation 

,, /drS(r.O) 
PS ~ JZ°dtfdrS(T,t) " 87R7 C ' K<6) 

if both conductivity and plasma frequency are known at a given temperature and magnetic field. 

In Bi-2212 crystals with T c w 85 K a c-axis conductivity a c w 0.5 ohm _1 cm _1 was found in Ref. 53 at T — 67pK 
and B z = 0.3 T. JPR frequency in a crystal with similar T c is uj p /2tt w 120 GHz at 4.2 K and w 90 GHz at 67 KB. 
From these data, using Eq. ( |5l| ) we estimate u) p /2n w 10 GHz at T = 67 K and B z = 0.3 T. With e c « 12 we obtain 
LOps/2-K Ri 0.2 GHz at B z = 0.3 T. Thus we estimate u> ps /ui p w 0.02 which justifies the static approach for calculation 
of plasma frequency. 

Consider now influence of the in-plane field on the c-axis conductivity. In the lowest order in Josephson coupling, 
we split [^„ )7l +i(0, 0) — Pn.n+i(Y, t)] into the contribution induced by pancakes and that caused by the unscreened 
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parallel component B x . Assuming that B x is along the x axis, we obtain a simple expression for the contribution of 
the parallel component to the phase difference: 

Vn,n+l(°>°) _ fn,n+l{r,t) w 

[<p n ,n+i(0, 0) - yn,n+i(r, t)] Bx =o - 2TTsB x y/<f> . (74) 
Inserting this expression into Eqs. (|7lj) and (|72|) we obtain 



(Tc(B z ,B x ) = (<K8J$/T) I drrG{r, B z )J (2nsB x r/<5> ), (75) 

where 

G(r,£y = / dtS(r,t,B z ). (76) 



The function S(r,t) describes the dynamics of the phase difference caused by mobile pancakes. If g(B z ,B x ) = 
Cc(Bj_, B x )/a c (B z ,0) is known in the vortex liquid, the correlation function G(r) — G(r) Jq $q/47tsT<7 c (.B z , 0) may be 
found using the inverse Fourier-Bessel transform: 

G{r,B z )= J dB x B x g(B x ,B z )3 (2irsB x r/^o)- (77) 
This procedure was realized in Ref. [53|and the function G(r, B z ) was found for pristine and irradiated Bi-2212 crystals. 



IV. LINE SHAPE OF PLASMA RESONANCE 

A. Formulation of the problem and general relations 

To find the line shape we need an equation for the phase difference which describes time variations of the interlayer 
phase difference and intralayer charges in the presence of an alternating electric field applied along the c axis. We 
will use the simplest version of this equation (p7|), which is similar to the Schrodinger equation. Thus the problem 
of inhomogeneous line broadening is very similar to the problem of the density of states of a quantum particle in a 
random potential, and one can use the same techniques. The role of the random potential in our problem is played 
by u n (r) — cos[(p^\ l+1 (r)} — C. Neglecting the Josephson coupling, we obtain: 

(u n (v)u m (0)) = ^-S nm . (78) 

where S (r) is the phase correlation function (||) . At large frequencies the random potential can be treated as a small 
perturbation. The exact criterion will be formulated below. Perturbative analysis of Eq. ( |27j ) can be conveniently 
performed using the Green function formalism. The Green function G mn (r,r';E) is defined as the solution of the 
equation: 

{ [E - C + u n {r)]5 nm + \ 2 jL nm V 2 } G mn . (r, r'; E) = S nn >8(r - r') . (79) 

I a 

with E — u> 2 /ujq. The homogeneous part of the alternating averaged phase is connected with G nn i (r,r';E) by the 
relation: 

where G(k, q; E) is the Fourier component of G nn i (r — r';E) = (G nn i (r,r';E)). The homogeneous part of the 
averaged electric field between layers n and n + 1 is: 

{E z {&)) - ^«„ +1 M) = ~- 2 G fo,0;S= 4) ^M- (81) 

Therefore, the dielectric function which describes the response to a homogeneous external electric field is connected 
with the Green function as: 

1 I^! G (Vo;£=4Y (82) 



e(uS) e c u\ 

The JPR absorption and shape of the JPR line is determined by Im[l/e(w)] 
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B. Diagrammatic expansion with respect to the random potential u„(r) 



The Green function at u(r) = is 

Go^i ^J^ , ^ q )=C + 2{i _J* 2 +s2/xlb (83) 

This function determines the spectrum of the phase collective mode w(k, q) = uj K(k, q). Using a standard diagram- 
matic approach, we can write: 

G(k, q; E) = [E- n 2 (k, q) - E(k, q; E)]" 1 , (84) 

where £(k, q;E) is the self-energy. In the lowest order with respect to w(r) (the Born approximation for scattering) 
we obtain: 

1 f AW da' 

£ (k, T,E) = -j -j^Go (k', q>; E) S(k - k'). (85) 

The imaginary part of the self-energy, S 2 (k, q; E), determines line broadening, while the real part, Si(k, q; E), deter- 
mines the shift of the resonance center due to inhomogeneities. For E — C <C A 2 /a 2 , one can neglect the k -dependence 
of 5(k) and we obtain: 

E 2 (k, q; E) = {S /16ir 2 ) J dkdqS [E-C- n 2 (k, q)] = (S /4\ 2 T )O{E - C), (86) 

where Q(x) = 1 at x > 0, and vanishes otherwise. So is the Fourier component of the function S(r) at k = 0, i.e. 

So = J dv (cos [^i 0) „ +1 (r) - <^ 0) „ +1 (0)] ) . (88) 

It is of the order a 2 in both the liquid phase and the vortex glass phase with strong disorder. In addition, in the 
liquid phase, the high temperature expansion gives the relation So = 2CT '/ 'Ej, which allows us to connect the width 
of the line with its center. Parameter E M = 5 /4A 2 7 = CT / (2E ) gives the reduced scattering rate due to the random 
potential. Perturbative expansion with respect to the random potential works if E — C ^> T.^. In the real units this 
criterion is given by (u 2 — ui 2 )/oj 2 3> T/2Eq. In spite of strong local variations of the random potential u(r), which are 
much larger than C, the parameter is much smaller than C. The reason is that typically the ac phase <fi' nn+ i(r) 
varies smoothly in space and effectively averages out these short-range variations. 
The dielectric function can be now represented as 

1 1 " 2 (»> 



Sl(M;£ ) = - " , n - (87) 



e(u) e c uj 2 - uj 2 [1 - si (u>)] - iuj 2 e(uj 2 - uj 2 ) ' 

Here = (T/2_Eo) 1//2 Wp is the line width for the vortex liquid at the right-hand side of the resonance line, and 
Si(ui) = (T/2ttE ) h\[co/a 2 (uj 2 — w 2 )]. The shape of the resonance line in the Born approximation is given by the 



p 

expression: 



Im 



1 



6( W ) 



^ - 2 e(c 2 -c 2 ), (90) 

£c(c 2 -c 2 [i- Sl H]) 2 + ^ 



b 



The Born approximation gives only the right-hand side of the resonance line because it describes transformation of 
the homogeneous plasma mode (with k = and q = 0) into the inhomogeneous phase collective modes with nonzero 
momenta k < a/ A 2 and q of the order unity in the presence of the random coupling cos^°| l+1 (r). Estimate for values 
k involved in the line broadening follows from the relation K 2 (k, q) w S , /4A 2 at values q of the order unity. This 
means that the main part of the JPR line near its center at u > uj p comes from averaging over distances R < A 2 /a 
in the vortex liquid state. Fluctuations of potential at smaller distances result in a far right-hand wing of the JPR 
line. Note that tot in the liquid phase is small compared with the frequency of the line center due to the small factor 
H = T/2E . 
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The absorption as function magnetic field is given by 



Im 



e(u,5) 



1 fiB r B 
T C [B- B r (l - Sl )f + faB r ) 



r 6 (B - Br) , 



(91) 



where B r — JujqEj&o/2Tu) 2 is the resonance magnetic field in the liquid phase. Eq. ( |9l| ) is valid at B—B r fiB r , and, 
thus the term (fiB r ) in the denominator may be dropped. Then the dependence of absorption on B is proportional 
to B j{B — B r ) 2 , and this dependence is in agreement with experiment! 2 ! 



C. Self-consistent Born Approximation 

Perturbative analysis of the previous Section can be improved by using the renormalized Green function (|84| ) in 
the Born expression for the self energy (|8q). This approximation, known as the self-consistent Born approximation, 
leads to the following equation for the self-energy 

Introducing dimensionless function si(£) and s 2 (C) °f the dimensionless variable £ = (E — C)/£oo + (2/tt) ln(A 2 /a 2 ) 
with Sqo = 5q/4Aj as: 



= s 2 (C) 5 v - = — ln ~r 



s 2 (C), ^=-z lll T2 L - s i(0 ( 93 ) 



we obtain the following equations for these functions 



S2 = — I — arctan (94) 

2 7T S 2 

. Sl = -i-ln[(C + S i) 2 +^] (95) 

The resonance absorption p(u>) can be written in terms of these functions as: 

p(o;K 2 /a; 2 = /[(a; 2 -^)/a; 6 2 ], (96) 

/(C) = [C + si(C)] 2 + ^(C)- (97) 

Here uj p is the resonance frequency shifted by the random Josephson coupling, a) 2 = oj 2 (1 — (T/ttE ) \ii(B/Bj)). 
This negative shift is due to the second order perturbative correction to the ground state "energy" , similar to the 
well known result of quantum mechanical perturbation theory. The scaling representation (pq) is exact, i.e., the 
shape of the line is fully determined by two dimensionless functions Si(£) and s 2 (C)- Eqs. (Q) and ( |95| ) allow us to 
calculate these functions approximately. Fig. ^| shows dependencies of the reduced JPR absorption p(u)w 2 / to 2 and 
the functions s\ t 2 on the reduced frequency £. The self-consistent Born approximation corresponds to summation of 
the perturbation series, which includes only certain class of diagrams. Omited terms have the same order and this 
approximation does not describe quantitatively the line shape in the region of maximum and below. In particular, 
within this approximation E 2 and resonant absorption still vanish at E < C — Yi^/n. In the real situation, there 
is a long absorption tail due to the localized plasma modes arising from exponentially rare fluctuation suppression 
of coupling at larg e areas. This tail is very similar to the Lifshitz tail in the density of states in inhomogeneous 
semiconductorsPir 5 !. 



D. Lifshitz Tail 



To calculate the exponential tail due,to localized plasma modes (left-hand side of the resonance line) one has to 
use the method of optimal fluctuationcJ. We will follow the derivations present in the book of Lifshitz, Pastur, and 
GredeskuEI. From now on we make the replacement E — > E — C. The resonant absorption is determined by the 
averaged spectral density 
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p(E) = i j dr' ^ Im [G„„, (r - r'; £ - tf)] 



which can be represented as: 

p(E) = J drJ2j Du(r)P{u n (r)}J2^A0,0;u)^ v (r,n;u)6(E-£ u {u}), 



(98) 



(99) 



where ^> v (r,n;u) and £ v are eigenfunction and eigenvalue of the state v in the system with the potential u n (r). 
The probability of realization of this potential is P{u n (r)} = exp[— H{u(r)}]. The scale of variations of typical 
eigenf unctions ^(r, n; u) for not too small £ u is much larger than the scale of variation of u n (r). This allows us to 
treat u n (r) as a short ranged Gaussian random variable and approximate H.{u n (r)} as 



H{u n (r)} = Y,J dr - 



(100) 



This is a good approximation to find the line shape not far away from the line center because the eigenvalue £q 
not far away from C comes from the effective random potential u n (r) averaged over large length scales R. This 
approximation breaks down only for very large and highly improbable fluctuations of the random potential. With 
exponential accuracy in the tail region we have p(E) cx exp(— &(E)) with 



$(£) = mmH{u n (r)}\ £o{uh 



(101) 



where £o{u} is the lowest eigenvalue at given u n (r). Following the standard line of reasoning^, which we outline in 
the Appendix, we obtain that $(£') is determined by solution of the nonlinear dimensionless eigenvalue problem: 



- V 2 ^„ + V 2 ^ = -e V^ n (r). 
with condition ipn{Q) = 1- $(E) can be found from solution of this equation as: 



HE) 



'04 



(102) 



(103) 



with 1P4 = J2 n I ^WnM- Numerical solution of Eq. ( |102| ) gives e = —0.1642 and ^4 = 2.4178. Asymptotic solution 
at large distances corresponds to a dipole field: 



V>n(r) = 0.224 



\e \r 2 - 2n 2 
(|e |r 2 +n 2 ) 5/2 



(104) 



because the time derivative of the phase difference <p n .n+i(r) = ^n(r) gives the electric field E z between layers n and 
n + 1. This electric field of the dipole type describes plasma oscillations with low frequencies when opposite charges 
are distributed between layers in the regions where Josephson coupling is below the average value due to fluctuations 
in the vortex den sity. 

From Eq. ( |l03| ) we obtain &(E) = 1A.72X 2 E/ So, which corresponds to the low-frequency tail (low-magnetic-field 
tail) of the dielectric function: 



Im 



1 




[<«>). 


cx exp 





-7.36 



E^wI-uj 2 ) 



To; 2 



(105) 



at fixed magnetic field (to p depends on B). At fixed frequency, the imaginary part of the inverse dielectric function 
as a function of magnetic field is: 



Im 



1 






oc exp ^ 



^B r — B 
pB r 



(106) 



where p = T/2E and the frequency dependent resonance field B r is defined after Eq. (^l|). In real units the coarse- 
grained optimum cosine fluctuation at E = C — u> 2 /u> 2 <C C is: 
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E I v 

u n (r) = cos^ ni „+i - C = - — rip n — 




Thus both sides of the resonance line are characterized by the width LOb « (T/ Eq) 1 / 2 
B z ^> Bj. Comparing the shape of the left-hand side of the resonance line, Eq. (1 
side, Eq. (BH), we see that the low- field side is-,sharper in agreement with experimen 



UJ, 



(107) 

in the vortex liquid phase at 
-with that of the right-hand 



By using Eq. ( |9l| ) and data of Matsuda et ala for the resonance line in the liquid phase at T — 36 K, we may estimate 
the parameter fi 0.06, while our theoretical estimate fi — T/2Eq(T) gives ~ 0.13 at X a b = 2500 A and s — 15.6 A. 
We note that, actually, Eq. (|9l]) may be invalid near B — B ri and, thus, our estimate of fi from experimental data is 
correct by order of magn itude only. The line was found to be anisotropic with a sharper left-hand side in agreement 
with Eqs. © and ( |l07| ). Thus, our model gives at least qualitative explanation of the JPR line form. 



V. CONCLUSIONS 



We have obtained a quite complete description of the JPR for the vortex liquid phase in magnetic fields B z ^> Bj. 
We related the JPR frequency in the presence of a tilted magnetic field, lo p (B x , B z ), to the correlation function of the 
pancake liquid. Numerical results for the density correlation function of the two-dimensional Coulomb plasma allow 
one to predict the field dependence of the plasma frequency or vice versa; experimental data for lo p (B Xi B z ) allow one 
to obtain information on the density correlations in the pancake liquid. 

We have explained the line shape of the JPR in the liquid phase by inhomogeneous broadening caused by disorder 
in pancake positions along the c axis. It is the excitation of the phase collective modes with nonzero momenta by a 
homogeneous-in-space alternating electric field in the presence of an inhomogeneous vortex distribution that leads to 
line broadening of the JPR line at a; > to p . Rare fluctuation suppressions of coupling in large areas due pancake density 
fluctuations result in the left-hand side (u> < u> p ) of the resonance line broadening. These mechanisms lead to a JPR 
line that agrees well with experimentally observed asymmetric line shape. Due to this inhomogeneous broadening of 
JPR line mechanisms associated with dissipation due to quasiparticles remain hidden when measurements are made 
in a disordered vortex state. 
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Appendix: Spectral density in the Lifshitz tail region 



In this Appendix we derive the eigenvalue equation (102) which determines the spectral density A(E) oc exp(— $(£')) 
in the Lifshitz tail region. Using the Lagrange technique we represent Q(E) as 



$(£?) = min [H{u n {r)} - /3{E - £ {u})} , 

u 

(3 is the Lagrange factor, which has to be determined by the condition £q{u} = E. Representing £o[u] as 



£o{u\ = min 



#o{*} + ]T J dvu n {v)^>l(v) 



with the Hamiltonian for the "free particle" 



(108) 



(109) 



(110) 



we can rewrite Eq. (101) as 



HE) 



-(3E + min /3H W} + min 

* \ u 



H{u n (r)} + [3J2 



dru n {r)^ 2 n (r) 



(111) 
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Optimal u n (r) is determined by the condition 



SH 



+ ^(r) = 



which gives 



and Eq. (101) can be represented as: 



Su n (r) 
u n (r) = -5 ^(r)/2, 

P 2 S Q 



$(E) = -0E + mm (/3H {^} - J dr <( r ) 



Optimal $„(r) is determined by the equation 

- A 2 7 ^L nm V 2 ^ m (r) - (/35 /2)*3(r) = £*„(r). 



(112) 



(113) 



(114) 



(115) 



Applying the operator — V 2 + s 2 /A 2 fc to this equation and neglecting small terms of the order s 2 /A 2 h we can rewrite 
this equation in a simpler form 



- A 2 V 2 *„(r) + (l3S /2)Vl*l(r) = -£V 2 „$„(r), 
where V 2 *„ = 2*„ - *„ + i - * n _i. Using relation 

H Q m = >iJ2l <kL nm VV n V* m = E / dr ((^o/2)<(r) + EV 2 n (r)) 



(116) 



(117) 



we obtain for Q(E): 



<Z>(E) = -PE + mm ( /3ff {*} ~ ^ E / dr *«( r ) 



/3 2 S 



E / dr *«( r )- 



Introducing dimensionless variables 



tf„(r) = # o (0)Vn(r), r = r 



2E 



v//35 v!/ 2 (0)/2' /95 *g(0)' 



we obtain the dimensionless equation ( 102 ). ^(i?) can be found from solution of this equation as: 

*(£?) = | A 2 * 2 (0)^ 4 



(118) 
(119) 

(120) 

(121) 



with ip4 = ^ n J dripn(r). Using the relation /3\l/g(0) = 2E/eoSo we obtain the final relation (103). 
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FIG. 1. Temperature evolution of pair distribution functions for a 2D vortex liquid (Langevin dynamics simulations). The 
plots are labeled by reduced temperature measured in units 2-kEq (as in Refs. pit). In parenthesis we also show approx- 
imate values of corresponding real temperature for Bi-2212 assuming a temperature dependent London penetration depth 
\ ab = 200nm/(l - (T/T c ) 2 ) 1/2 with T c = 90K. 
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FIG. 2. Phase correlation functions at different temperatures obtained by numerical integration of Eq. (^) with numerical 
pair distribution functions shown in Fig. |l|. 
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FIG. 3. The Fourrier transform of the phase correlation function f 3 (b) which determines the angular dependence of the 
plasma frequency given by Eq. (R). The function is calculated by the Fourier-Bessel transformation (|52j) using the phase 
correlation functions shown in FigTM 
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FIG. 4. Comparison between the function f(b) = 
(optimally doped Bi-2212 sample at 40 K from Ref. 



(6)// s (0) extracted from the angular dependence of the resonance field 
Fig. 3) and numerically calculated f(b) at T = 0.032 • 2nE 
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FIG. 5. Temperature dependence of the phase correlation length defined by Eq. ( ff6| ) in units of lattice constant ao. Squares 
represent results of numerical integration of Eq. (^) using pair distribution functions h(r) generated from Langevin dynamics 
simulations. Spread of data is due to an error in the integration because of numerical noise in h(r) at large r. Two filled 
circles represent correlation lengths extracted directly from the phase correlation function, which was obtained by Monte Carlo 
simulations of the frustrated XY modelE3. 





FIG. 6. Plot of reduced JPR lineshape p(ui)u b /uj vs reduced frequency £ obtained within SCBA (Eqs. 
The insert shows frequency dependencies of the dimensionless self energies si and S2. 
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